Interface roughening with nonlinear surface tension 
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Using stability arguments, this Brief Report suggests that a term that enhances the surface 
tension in the presence of large height fluctuations should be included in the Kardar-Parisi-Zhang 
equation. A one-loop renormalization group analysis then shows for interface dimensions larger than 
~ 3.3 an unstable strong-coupling fixed point that enters the system from infinity. The relevance of 
these results to the roughening transition is discussed. 
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During recent years, kinetic roughening of growing in- 
terfaces has been the object of intense research. In 1986, 
Kardar, Parisi, and Zhang (KPZ) jl| suggested the fol- 
lowing Langevin equation to describe the macroscopic 
dynamics of a stochastically growing interface: 



(1) 



Here, h = h(x,t) is the (coarse-grained) height of the 
d-dimensional surface, and r\ is a stochastic noise that 
roughens the interface. It is uncorrelated and Gaussian 
distributed, i.e. {rj(x,t)) = and 



{r)(£, t)r,{x',t')) = 2DS d {x - x')5{t - t') 



(2) 



The first term on the right-hand side of Eq. ([!]) is a sur- 
face tension that tends to smooth the interface; the sec- 
ond term accounts for growth perpendicular to the sur- 
face orientation. In principle, a constant term has to be 
added, which, however, can be absorbed in the definition 
of the height variable. 

A comprehensive discussion of Eq. ([!]) can be found 
in Near the stationary state, the height profile 

is self-affine, i.e. invariant under a rescaling h'{x,t) = 
b~ x h(bx, b z t). The exponent \ is the roughness exponent, 
z is the dynamical critical exponent. Eq. ([!]) is invariant 
under the transformation h'(x, t) = h(x + Xvt, t) + v ■ x, 
that tilts the interface by a small angle v (J], leading to 
the exponent relation 



x- 



z = 2. 



(3) 



In one dimension, the critical exponents \ an( i z are ex ~ 
actly known, since the stationary solution of the Fokker- 
Planck equation 



dP 
~di~ 



is given by 

P[h(x)} = exp 
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leading to \ — 1/2 and (with Eq. (||)) z = 3/2 (see e.g. 
I- 

In higher dimensions, the stationary distribution is 
not exactly known, and the analytical approaches that 
have been taken so far, are renormalization group the- 
ory illfl@] an d self-consistent methods A renor- 
malization group calculation to one-loop order gives the 
following flow equations for the parameters J3J" 



Ml) m f 2-d 
dD{l) 



dl 
d\{l) 

dl 



D(l)(z-d-2 x + g) 
X(l)( x + z-2), 



(5) 



where g = K^A A D/4v 3 is an effective coupling 
constant, with Kd being the surface area of the d- 
dimensional unit sphere, divided by (27r) d , and A the 
cutoff for the wave vector. Equations (o) can be com- 
bined to form an equation for g alone, 



dg 
dl 



, 4d-6 
d+^^g 



(G) 



In one dimension, this equation has a stable fixed point 
at g* = 1/2. Inserting the fixed-point value in the flow 
equations for v and D gives again the above-mentioned 
values for the critical exponents. The second fixed point 
g* = is unstable. 

In 2 dimensions, the fixed point g* = is still unstable, 
and there is no stable fixed point at all. In dimensions 
larger than 2, the fixed point g* — becomes stable, and 
an unstable fixed point g* = d(d — 2)/(4d — 6) appears 
that separates the region where the flow goes to zero from 
the region where it diverges. These results indicate that 
there exists a phase transition from a flat interface to a 
rough interface for dimensions larger than 2. Although 
the rough phase is presumably controlled by an attrac- 
tive fixed point, this fixed point is not accessible within 
a one-loop renormalization group. A two-loop calcula- 
tion yields similar results 0, and recently it has been 
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shown that a finite stable fixed point is not reached by 
perturbative calculations to any order |l4| . 

There are certainly many possible explanations for this 
failure of perturbation theory (e.g. lack of a small param- 
eter, or nonanalytical behavior). In this Brief Report, I 
suggest a further explanation by noting that Eqs. (|5|) 
might reflect a physical instability of Eq. ([[]) , and that a 
stabilizing term is necessary. Integrating the flow equa- 
tion for v (Eq. (||)) in the neighborhood of v = 0, one 
finds that v flows through zero and becomes negative for 
some finite value of / (see also pi ). If this behavior 
reflects the physics of the system, it indicates that the 
interface is unstable, since the surface tension becomes 
negative. A self-consistent mean-field theory Jl6| for the 
KPZ equation is also unstable, suggesting an instability 
at least in high dimensions. A stable mean-field theory 
is presented in ITj]], where a nonlinear surface tension 
is added to Eq. ([l]). This mean- field theory leads to an 
interface that has bumps of a characteristic height. Insta- 
bilities have even been found in computer simulations of 
the (discretized) KPZ equation in two dimensions [fl7|| . 
A continuum limit, however, is only well defined if the 
discrete model leads to a stable and smooth interface at 
large scales, and it seems therefore necessary to include 
a stabilizing term in Eq. ([!]). 

The equation studied in this paper is the following [O : 



9h „ 2l 



T). 



(7) 



The new term proportional to k is generated at the sec- 
ond order in the expansion for the surface tension (see e.g. 
|^|) and is the stabilizing term in the above-mentioned 
mean- field theory |T?{] . It is possible that the nonlinear 
surface tension, although present at microscopic scales, 
becomes irrelevant at larger scales, which is indeed the 
case in the neighborhood of the fixed points g* = and 
g* = d(d — 2) /(Ad — 6), as we shall see below. We will, 
however, find an additional fixed point where k is not 
small, suggesting that the stabilizing term is important 
even at large scales in the rough phase. 

Let us first discuss the effect of the new term on the 
one-dimensional KPZ equation. The stationary height 
distribution is now 
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Under a rescaling x — > x' = x/b and h — » hi = h/b x with 
X = 1/2, the second term is multiplied by 1/6, while the 
first term remains invariant, indicating the irrelevance 
of k/D. A one-loop renormalization group calculation 
for the stationary distribution in Eq. (||) leads to the 
same conclusion. In momentum space, the propagator is 
D/vk 2 , and the vertex is — (n/2AD)k 1 k2k' i k i 5(ki + k2 + 



&3 + ki). To one-loop order, one obtains after a short 
calculation the following flow equations: 



d{u/D) _ v 
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leading to 



dc 
1l 



c(-l - 5c) 



(9) 



(10) 



for the effective coupling constant c = (kD /2v 2 )KiA. 
The fixed point c* — is stable, indicating that the stabi- 
lizing term can be neglected at large length scales. There 
is also an unstable fixed point at c* = —1/5. For initial 
values c < —1/5, the flow of c goes to — oo, indicating 
that the surface becomes unstable in this parameter re- 
gion. We will not consider further the region of negative 
c. Since the negative fixed point must also occur in the 
flow equations found from dynamics, Eq. (^), it can how- 
ever be used to make sure that the calculations do not 
contain errors. 

Let us now renormalize the equation of motion Eq. (|t]) 
to one-loop order. Inserting the equation of motion in 
the Gaussian probability distribution of the noise 



W[rj] oc exp 



1 

Id 



d d x / dtri(x,tf 
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and introducing an auxiliary field h, we obtain the weight 
of a given space-time configuration [h(x, t)] Jig] 

W[h] oc J V[ih]exp{j[h,h]} , 

with the dynamical functional 



J[h, h] 



d d x 



Dhh — h x 



^ _ v ^ h - ^(v^) 2 - -V 2 h(Vh) 2 
dt 2 K ' 2 x ' 



(12) 



The dynamical functional J plays the same role in dy- 
namical renormalization group as the Hamiltonian in 
statics. The propagators of this model are 



G (k, t) = {h(~k, t)h(k, t)) = 6{t)e 



-uk A t 



and 



Cb(M) = (h(-k,t)h(k,t)) = De- vk2 ^/vk 2 
and the vertices are 

_ ^(^2 • hjh^h^h^Siki + k 2 + k 3 ) 

and 



(13) 
(14) 
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-k 2 (k 3 • ^)Hh%Ji% 3 h%J{ki +k 2 + k 3 + fc 4 ) . 

Here, (. . .}o indicates an average within the Gaussian the- 
ory, where the nonlinear terms proportional to A and 
k are not taken into account. Renormalization of this 
model is done by first integrating over the large wave 
vectors A/b < k < A, where A is the wave vector cutoff, 
and b = 1 + 1 is close to 1. Next, the system is rescaled 
to the original size by introducing new variables k' = bk, 
t' = t/b z , h' = h/b x . To one-loop order, this gives the 
following flow equations for the parameters: 



dv ( X 2 D 2 - d T , , d 9 kD „ k , 

-dl=»{ Z - 2 + ^— K « A +2^ A 
dD ( \ 2 D , . , 



d=2 



dX , / „ k-D T „ 4 rf 



k(^2x + z-4- — — — K d A a 



dl 



X 2 D 6 - d 
Av 3 d 



2v 2 d 
K d A^ 2 



(15) 



and finally for the effective coupling constants g = 
X 2 DK d A d ~ 2 /Av 3 and c = nDK d A d /2v 2 , 



dg ,„ , 4d-6 

dc , , 4d - 10 d + 4 . 
dl =c{ - d+ ^^ g -— c) 



(16) 



The diagrams contributing to k that contain four A 
vertices, cancel. This is a consequence of the tilt in- 
variance of the KPZ equation: If k = 0, the system is 
tilt invariant, and this property must not change under 
rescaling. Therefore the flow diagram on the <?-axis looks 
exactly as in the absence of k. There is consequently 
a fixed point at g* = c* = that is unstable in one di- 
mension, and stable above d = 2. The fixed point c* = 0, 
g* = d(d~ 2)/(4d— 6) is stable in d < 1.5, and is a saddle 
point for d > 1.5. A further fixed point is given by the 
intersection of the lines 2 — d + g(4d — 6)/d — 5c = and 
-d + g{Ad - 10) /d - c(d + 4) /d = 0. For d > 3.295, this 
fixed point lies in the physically interesting region g > 0, 
c> 0, and is given by g* = d(2d 2 -d+4)/(8d 2 -30d+12) 
and c* = (6d 2 - 10d)/(8d 2 - 30d+ 12). A linear stabil- 
ity analysis reveals that this fixed point is unstable, with 
two complex eigenvalues in dimensions between 3.574 and 
11.53, and with two positive eigenvalues in all other di- 
mensions. For dimensions close to 3.295, this fixed point 
is very large, indicating that it does not split from one 
of the other fixed points, but that it enters the system 
from infinity. The qualitative flow diagrams in one, two, 
three, and four dimensions are shown in Fig. |l|. 




g g 

FIG. 1. Schematic flow diagram obtained from Eq. (16) in 

one, two, three, and four dimensions. 



The values of the critical exponents z and x a t the new 
fixed point are z = (2d 3 + 5d 2 - 44d+ 16)/(8d 2 - 30d+ 12) 
andx= (-2d 3 + 17d 2 -26d + 8)/(8d 2 -30d+12). This 
one-loop result is certainly still far from reflecting the 
true behavior of the system, with a stable strong-coupling 
fixed point that occurs already in two dimensions, and 
with < x < 1 ancl z < d. Nevertheless, we can draw 
some interesting conclusions: (i) In the presence of the 
stabilizing term, a fixed point can be found in the strong- 
coupling regime. Since this fixed point is fully repulsive 
in one-loop approximation, its presence does not affect 
the behavior of the system at large scales. However, it 
is possible that this fixed point becomes attractive and 
occurs already in two dimensions, when a better approx- 
imation is used, (ii) The strong-coupling fixed point en- 
ters the system from infinity and is not close to the g axis. 
If this is not just an artifact of the perturbation theory 
but reflects the true behavior of the system, there is no 
reason to expect that other higher-order terms can be ne- 
glected, (iii) The stabilizing term and other terms that 
do not preserve the tilt invariance of the KPZ equation, 
cannot be generated by the renormalization procedure. 
They have to be included explicitely in the beginning, 
(iv) The mapping of the KPZ equation on directed poly- 
mers in random media via the Cole-Hopf transformation 
cannot be performed any more in the presence of 
the nonlinear surface tension. Nevertheless, more direct 
mappings between growing Eden clusters and directed 
polymers |l!| suggest that the two systems should be 
equivalent even when higher-order terms are taken into 
account. Alternatively, as proposed in |l7j], one might 
ask the question if the strong-coupling behavior can at 
all be characterized by a single universality class. 

Calculations to higher loop orders or the inclusion of 
further terms might produce better results than the one- 
loop calculation presented here, but they are uncontrolled 
as long as the size of the neglected terms cannot be esti- 
mated. Self-consistent methods might be more success- 
ful than perturbation theory. Mode-coupling theory has 
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been shown to yield excellent results in one dimension, 
since the neglected vertex corrections are small |l^,[l3| . It 
still has to be seen if this holds also in higher dimensions. 
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